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Abstract 

In the mathematical problem of linear hydrodynamic stability for 
shear flows against Tollmien-Schlichting perturbations, the continuity 
equation for the perturbation of the velocity is replaced by a Poisson 
equation for the pressure perturbation. The resulting eigenvalue prob- 
lem, an alternative form for the two - point eigenvalue problem for the 
Orr - Sommcrfcld equation, is formulated in a variational form and 
this one is approximated by finite element method (FEM). Possible 
applications to concrete cases are revealed. 
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1 The classical model and the model based on a 
Poisson equation for the pressure perturbation 

The mathematical problem of linear hydrodynamic stability is [1] 
du' 

— + (U • grad)u' + (u' • grad)U = (la) 

= -grad p 1 + Re' 1 Au', (t,x) el + x!l, 

div u = 0, (i,x) £l + x!l, (lb) 

u'(t,x) = 0, t G R+, x G dfl, (lc) 

u'(0,x) = u' (x), x G fi, (Id) 

and its solution is the perturbation (u',p f ), for t > 0. Here all variables are 
non dimensional, (U, P) is the basic motion in the domain £1 C R n , n = 2 
or 3, Re = U^L/v is the Reynolds number, L is a length scale, Uoo is a 
velocity scale, v is the coefficient of kinematic viscosity. 
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Applying the divergence operator to (|lap and using the divergence-free 
condition on u', we obtain 



Ap' 



div ■ ((U • grad)u') — div • ((u' • grad)U) 

(t,x) g r + x n, 



(2) 



and projecting (fla|) on n*, we have 



n*(x) • gradp'(t,x) = n*(x) • F(Re, U(x), u'(t, x)) 
(£K + , x G <9fi, 



(3) 



where n* is the unit outward normal to <90 and F(i?e, U, u') = i?e _1 Au' 
— (U • grad)u' — (u ; • grad)U. Another form of the right hand - side of 
the Poisson equation for p' can be obtained by using the divergence-free 
condition on U. 

Thus the mathematical problem of the linear hydrodynamic stability 
becomes (fTa|) . ([2]), (fTcj) . Q, (|ld|) . Of course, it is assumed that equation 
(fTa|) can be continued on the dQ. 

2 An alternative approximation of the two - point 
eigenvalue problem for Orr - Sommerfeld equa- 



Let Q, = {x = (x,y,z) G R 3 \ - oo < x,z < oo,0 < y < a} (a > 1) 
and assume that the basic flow is of the form U(x, y, z) = (U(y), 0,0). 
Let us choose Tollmien-Schlichting waves - like perturbations u'(t, x, y, z) = 
u' (y)exp[ia(x-ct)}, u' (y) = (u(y),v(y),0), u'(t,x,y,z) = (u'(t, x, y), v'(t, x, y), 0), 



p'(t, x, y, z) = p(y)exp[ia(x — ct)], where i = a is the streamwise wave 

number, c = c r + the wave speed and q is the amplification rate. 

In this case, model ([T]) leads to the classical two - point eigenvalue prob- 
lem for Orr - Sommerfeld equation in (c,<p), where <p is the nonexponential 
factor of the stream function, while (fTa|) . ([2]), (fTc|) . ([3]), (jldp becomes the 
two - point eigenvalue problem in (c,u,v,p) 



—aRe 1 iu + (aRe) 1 iu" + Uu — a 1 iU'v+p = cu, y G (0, a) ,(4a) 



where the prime stands for the differentiation with respect to y, c is an 
eigenvalue and (u,v,p) is an eigenvector. 



tion 




—aRe iv + (aRe) 1 iv" + Uv — a 1 ip' = cv, y G (0, a) 

—a 2 p + p" = —2aiU'v, y G (0, a), 

u(y) =0 for y = 0,y = a, 

v(y) =0 for y = 0,y = a, 

p'{y) = Re~ l v"{y) for y = 0, y = a. 



(4b) 
(4c) 
(4d) 
(4e) 
(4f) 
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Let L 2 (0, a) be the Hilbert space of all measurable complex-valued func- 

a 

tions u, defined on (0,a), for which J \u(x)\ 2 dx < oo and the inner prod- 

o 

a 

uct is (u, v) = J uvdy. Denote Du the generalized derivative and consider 
o 

the spaces H%(0,a) = {u G L 2 (0,a)\Du G L 2 (0,a),u(0) = u(a) = 0}, 

a 

Lg(0,a) = {p G L 2 (0,a) | / p(x) da; = 0}. 

o _ _ 

Multiply (|la|) . (|4b|) and (ficj) by arbitrary test functions / 1; / 2 and ~g 
respectively, integrate over (0,o), apply partial integration if necessary and 
take into account the two point conditions (|4dj) - (j4f|) to obtain the following 
weak formulation of problem in (c,u,v,p) G C x Hq(0,o) x Lq(0, a) 

a a a 

— aRe~ x i j uf^dy — (ai?e) _1 i ^ u'/j dy + j Uufidy— 



o oo 

a a a 



(5) 



a~ l i j U'vf 1 dy + Jpf 1 dy = c- J ufidy,V fi G Hq(Q,o) 



a a a 

-aRe~ x i J v f 2 dy - (aRe)" 1 i J v'J 2 'dy + J Uvf 2 dy 
o oo 

a a 

+ a _1 i J pf 2 'dy = c- J vf 2 dy,V f 2 G H^(0,a), 



(6) 



(7) 



a a 

- a 2 J Pgdy - J p'g'dy = -Re~ l v"(a)g{a) 
o o 

a 

+ Re~ 1 v" (0)5(0) - 2m ^ U'vgdy,V g G Lo(°> a )- 



3 Approximation of problem fl5]) - (|7j) by FEM 

In order to perform this approximation, let us divide the interval [0, a] 
in N + 1 subintervals K = Kj = [yj, yj+i], < j < N, where = yo < Vl < 
• • • < Un+1 = The sets K represent a triangulation Th of [0, a]. 

The approximate basic shear flow \Jh(x,y,z) = (Uh(y), 0, 0) is defined 
by Uh(y) = U(y), y G [0, a]. The approximate amplitudes Uh(y), Vh(y) and 
correspond to the exact ones, v(y) and p(y), respectively. 
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Introduce the spaces Vh = {v : [0, a] — > C | v G C[0, a], v is on i£j 
an one - dimensional polynomial, having complex coefficients, of degree 2, 
< j < N, u(0) = u(a) = 0} and M h = {v : [0, a] -> C | u G C[0, a], v is on 

an one - dimensional polynomial, having complex coefficients, of degree 
1, < j < N}. 

Correspondingly, variational problem ([5|) - (J7|) is approximated by the 
following problem in (c h ,u h ,v h ,p h ) G C x V h 2 x M h , 



a a a 

aRe~ l i J u h f x dy - (aRe)~ 1 i J u' h J 1 'dy + J XJu h J x dy- 
o oo 

a a a 

- a-h J U' Vh J ld y + J Ph - h dy = c h • / u^dy, V h G V h 







ai?e l i j v h f 2 dy - (a Re) 1 i J v' h f 2 dy + / Uv h f 2 dy 
o oo 

a a 

+ a~ l i J pij 2 dy = c h ■ j v h f 2 dy,V f 2 G V h , 



(9) 



.2 



a 
o 



y P/>5^ - J P h g'dy = -Re 1 v' h \a)g(a) 
o o 

+ i?e- 1 <(0)5(0) - 2m / U'v h gdy,V g G Af ft . 



(10) 



In order to obtain u^, v/i, we use a basis of real functions of Vh- Let 
J = Jk = {1,2,3} be the local numeration for the nodes of K, where 1,3 
correspond to yj,yj + \ respectively and 2 corresponds to a node between yj 
and yj+i- Let {(/>„, n G J} be the local quadratic basis of functions on K 
corresponding to the local nodes. Let J* = {1, . . . , 2N + 1} be the global 
numeration for the nodes of [0, a], where the nodes corresponding to yo and 
yN+i axe not taken into account, and let L\ be a matrix whose elements 
are the elements of J*. Its rows are indexed by the elements K G Th 
and its columns, by the local numeration n G J. We take the value at 
the locations of L\ which we do not consider in the computations, i.e. the 
locations where K = Kq, n = 1 and K = Kn, n = 3. Write ra* = L\{K,n) 
or, simply, n t for the element re* of L\ which depends on K and n. Let 
{(^n*, 0), (0, $n„); n* G J*} be a basis of functions of V^. If n* corresponds 
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to yj, then $ n<< is a real quadratic function on Kj—i and Kj and its value 
is zero on [0, a]\(Kj-i U Kj). If n* lies between yj and j/j+i, then $ nt 
is a real quadratic function on Kj and its value is zero on [0, a]\Kj. We 
have = (j> n (y), where y £ K, = L\{K,n). Let u nt , v nt be the 

values of Uh, Vh at the nodes n*, n* G J*. Retaining our convention about 
= Li(K,n), we do not write in the sequel the conditions n ^ 1 for 
X = i<Lo and ?i / 3 for .fT = .Kjv- We have 

n*gj* i^G T h ,yeKn& J 

and a similar expression for Vh(y). 

In order to obtain p^, we use a basis of real functions of M^. Let 
/ = Ik = {1, 2} be the local numeration for the nodes of K, where 1, 2 cor- 
respond to yj,yj + i respectively. Let {ip m ,m £ /} be the local affine basis of 
functions on K corresponding to the local nodes. Let I* = {0, 1, . . . , N, N + 
1} be the global numeration for the nodes of [0, a] and let L2 be a matrix 
whose elements are the elements of I*. Its rows are indexed by the ele- 
ments K G 7ft and its columns, by the local numeration m £ I. Write 
m* = L2(K,m) or, simply, m* for the element m* of L2 which depends 
on K and m. Let {^m^m* G /*} be a basis of functions of M^. If m* 
corresponds to yj, then ^> m , is a real affine function on Kj_\ and and its 
value is zero on [0, a]\(Kj-i UKj). We have * m ,(y) = ^ m (y), where y G K, 
m* = L2(K,m). Let p mt be the values of p/j at the nodes to*, to* G /*. We 
have 

Ph(y) = Prn^m*(y)= ^ Yl PmAm{v) ■ 

m*Gi» Ke T h ,y£Km£ I 



Problem (J8J) - ()10p becomes the following algebraic eigenvalue problem 
in (c h , «i, vi, ... , U2JV+1, W2JV+1, Po, • • • , PJV+i) G <C 1+2 ( 2Ar+1 ) +JV + 2 



u nt [—aRe l i j (pn^^dz — (aRe) l i j (f)' n (f>' k dz+ 
Ke T h n£ J K K 

+ j U(f) n (pkdz]+ ^2 [-a^i J U'(p n (f) k dz]+ ^ 

K A'G T h n£ J K 

+ ^2 Y Vm * j ^^kdz = C h -{ ^ Y Un * J <Pn<f>kdz} 



Ke Th,me I K Ke T h ne J K 
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j k dz 



(13) 



^ ^ Un.f-aiie x i y 4> n <f) k dz - (aRe) 1 i j <j)' n (p k 

Ke T h n£ J R K 

K Ke Th me I K 

= c h ■ { X] Y Vn * \ ^kdz} 
Ke T h ne J J K 

Y Pm * [ _Q;2 / ^m^edz - J ip' m ip' e dz] 
Ke T h me I K K 

= - v n ,Re- l ^{a)^ N+1 {a) 

K=K N ; ne Jpc N ;n^3 

+ Yl wn,fle -1 <.(0)*i(0) 

K=Ko; ne JK ,nj^l 

+ Y Y Vn *(~ 2ai / U'(/) n ipedz), 
Ke Thne J J K 

for all k £ J, £ £ I, for all K G Th, k ^ 1 when A" = K , 7^ 3 when 
K = K N . 

The eigenvalue problem in G C, = ( ifi,«i, . . . , U2AT+ij V2N+1) £ 
C 2(2iv+l) ; = ( p , Plj . . . , PN+1 ) £ C N+2 represented by relations §XB) - 
()13p can be written in the following matrix form 

K h A h + L h B h = c h S h A h , G h B h = H h A h , (14) 

where K h ,S h G C 2 ^ 2N+1 ^ 2N+1 \ L h G C^ 2N+1 ^ N+2 \ G h G C ( N+2 ^ N + 2 \ 

H h £ £{N+2)*2{2N+l) _ 

Matrix is positive definite and symmetric. So problem (|14p becomes 

(i^ + L h Gl x R h )A h = c h S h A h , B h = G^H h A h , (15) 

The eigenvalue problem given by the first equation of (|15p can be solved by 
using the QZ method or, after multiplication by Sj^ 1 at the left, by using 
the QR or LR method. 

4 Possible applications to concrete fluid flows 

Once (|15p solved, several linear hydrodynamic stability characteristics 
can be determined numerically. Keeping a and Re fixed, (|15|) yields ap- 
proximations of c and of amplitude distributions u(y), v(y) and p(y). Re- 
peating this operation with various values of a and Re, we can identify the 
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pairs (Re, a) where the approximate value q = 0. Thus, we can construct 
the neutral curve. We also can determine the curves Ci(Re,a) = constant 
of constant amplification factor and the wave speed c r (Re, a), for (Re, a) 
belonging to the neutral curve. 

In particular, problem (|15p is the basis of a computer program for the 
case of Prandtl's boundary layer, which will be presented elsewhere. 
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